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Abstract. One of the so-called viable modified gravities is analyzed. This kind of gravity theories are characterized by a well 
behavior at local scales, where General Relativity is recovered, while the modified terms become important at the cosmological 
level, where the late-time accelerating era is reproduced, and even the inflationary phase. In the present work, the future 
cosmological evolution for one of these models is studied. A transition to the phantom phase is observed. Furthermore, the 
scalar-tensor equivalence of f(R) gravity is also considered, which provides important information concerning this kind of 
models. 
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INTRODUCTION 

As a consequence of the discovery of the accelerating ex- 
pansion of the universe, a lot of dark energy candidates 
have been proposed along recent years. Modified grav- 
ity theories has become one of the most popular candi- 
date, since it provides a way to understand the problem 
of the dark energy and the possibilities to reconstruct a 
gravitational theory capable to reproduce late time ac- 
celeration, and even the inflationary phase (for a review 
see Ref. [1]). In this sense, f(R) theories are the simplest 
ones, since they consist in a generalization of the Hilbert- 
Einstein action to a more complex function of the Ricci 
scalar, 



S = jd 4 x^j[f(R)+2K 2 J? m ] 



(1) 



Here the coupling constant is as usual given by K — 
8kG, while Jzf,„ is the matter Lagrangian. The field equa- 
tions corresponding to the action (1) are obtained by the 
variation of this action with respect to the metric tensor 
gjiv, what yields 
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Here the subscript r denotes derivatives with respect 
to R. Here we are interested to study flat Friedmann- 
Lemaitre-Robertson- Walker (FLRW) metrics, 



ds 2 = 
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Hence, modified FLRW equations are obtained through 
the field equations (2), 
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where dots denote derivatives with respect to the cos- 
mic time, H(t) = a/a is the Hubble parameter, and Ricci 
scalar is given by R = 6 (2H 2 + H) for the metric (3). 
Hence, any cosmology can be reproduced by the ap- 
propriate action (see Ref. [2]). In this manuscript, we 
focus on the study of a subclass of modified gravities, 
the so-called viable f(R) gravities, as they accomplish 
some indispensable conditions to be considered realis- 
tic candidates to describe the universe evolution. This 
kind of f(R) gravities are usually described by actions 
of the type f(R) = R + F(R), which basically repre- 
sents Hilbert-Einstein action plus an additional term that 
should affect only at cosmological scales, while at lo- 
cal scales, General Relativity is recovered. Viable grav- 
ities are able to satisfy this constraint, as well as avoids 
large instabilities in the presence of matter distributions, 
and the anti-gravity regime, by the appropriate form of 
the function F(R) (see [3]-[4]). Moreover, these theories 
are capable to reproduce the effects of dark energy, and 
even inflation. Note that both equations in (4) are writ- 
ten in such a way that F(R)-teims are on the matter side. 
We can define an effective energy density for the extra 
F (R)-terms, so that the first FLRW equation in (4) can 
be rewritten for the kind of actions in terms of the cos- 
mological parameters £1,-, 



i — n„, + £if(R) , 



where, 
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Here, we have assumed f(R) = R + F(R). Then, the first 
Friedmann equation (4) takes a simple form, with two 
fluids contributing to the scale factor dynamics. In ad- 
dition, the continuity equation V^T^ = for a per- 
fect fluid with an equation of state (EoS) p m = w m p m 
yields p m + 3H(1 + w m )p m = 0. Hence, for a particular 
F(R), the corresponding cosmological evolution can be 
obtained through equations(4), and the continuity equa- 
tion. Then, under these circumstances, the F(R) term 
contributes to the dynamics of the scale factor, and it is 
capable to reproduce both accelerating epochs of the uni- 
verse evolution. Here we will focus on the study of one 
model of this kind of viable modified gravities, proposed 
by Hu and Sawicki in Ref. [3], and whose action is given 

by, 

F H s(R)—RHS c2{R/RHs)n + l , (6) 

where {c\,C2,n} are free parameters and Rhs = Kl P m - 
The next sections are devoted to the study of the cos- 
mological evolution for this model, specially during the 
dark energy epoch, and in the future. We will also study 
the behavior of their scalar-tensor counterpart, where the 
phase space is explored. This analysis is based on the 
study performed in Ref. [5]. 



COSMOLOGICAL EVOLUTION 




FIGURE 1. Evolution of the EoS parameter (7) as a function 
of the redshift z for the model (6) with n = 1, c\ = 2 , C2 = 1, 
and initial conditions h(0) = 1 and h'(0) = 0.4. 




In this section, we explore the cosmological evolution 
for the Hu-Sawicki model. For convenience, we express 
the equations (4) in terms of the redshift z instead of the 
cosmic time f, 1 + z = ^j, where «o is the value of the 
scale factor at the present time to, such that the current 
epoch corresponds to z = 0. 

Then, we can fit the current values of the cosmolog- 
ical parameters using the observational data [6], where 
H = 100/z km s- 1 Mpc- 1 with h = 0.71 ±0.03 and 
the matter density, QP m = 0.27 ± 0.04, while the matter 
fluid is considered pressureless (cold dark matter and 
baryons), w m = 0. Then, by setting the initial condi- 
tions, the corresponding cosmological evolution can be 
obtained through the Hubble parameter in terms of the 
redshift and the future evolution can be explored (with 
— 1 < z < 0). Here, we analyze the evolution of the effec- 
tive EoS w e ff, which is defined as 
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FIGURE 2. Evolution of the EoS in the interval - 1 < z < 
—0.8, plotted in more detail, where the EoS crosses clearly the 
phantom barrier and oscillates. 



For simplicity we redefine the Hubble parameter as 
H(z) = Hq h(z), such that h(0) = 1, while the initial con- 
dition on the first derivative h'(0) can be fixed by assum- 
ing that our model mimics ACDM at z = 0, which basi- 
= |£2„. Let's now analyze the 



(7) 
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Fhs(R) model in (6), assuming c\ = 2 and C2 = 1, which 
are dimensionless parameters, and a power of n = 1 . In 
Fig. 2, the evolution of the EoS parameter (7) is plotted 
with respect to the redshift for the initial conditions that 
mimic ACDM model at z = 0. In spite of the evolution of 
the EoS parameter is quite similar to the ACDM model 
for positive redshifts, the universe enters in a phantom 
phase in the future (negative redshifts), where the EoS 
parameter turns out less than — 1, and it oscillates close 
to z = — 1 as shown in more detail in Fig. 2. In Fig. 3, 




0.0 0.5 



FIGURE 3. Evolution of the cosmological parameters 
{£l m ,Slp} defined in (5) as a function of the redshift z for the 
Hu-Sawicki model. 



the evolution of the cosmological parameters {£2„,,n/r} 
with respect to the redshift is plotted, where Q.f tends to 
dominate completely the universe in the future. Note that 
the viable modified model considered here (6) produces 
some oscillations along the cosmological evolution. In 
general, the transition to the phantom epoch occurs for 
this class of models, as pointed out in [5], which does not 
imply the existence of a singularity in the future. Never- 
theless, a singularity may occur along the whole evolu- 
tion, since the scalar potential presents a pole for a finite 
value of the scalar field. However, the singularity may be 
excluded from the range of interest by a particular choice 
of the free parameters. 



SCALAR-TENSOR COUNTERPART 



where we have assumed n = 1 in (6), as in the previous 
section. Then, the scalar potential (9) yields, 
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Hence, the scalar-tensor representation is not uniquely 
defined, but the potentials exhibit two branches, which 
in principle do not affect the cosmological evolution (see 
Fig. 4), but it may influence the behavior of the phase 
space. In addition, the model introduces a boundary con- 
dition on the value of the scalar field, being <p < 1, a 
limit where both branches of the potentials converge, and 
where the first derivative of the potential exhibits a pole, 
which induces a sudden singularity (see Ref. [8]). Note 
that there is a direct correlation between the behavior of 
the scalar field and the evolution of the EoS parameter 
as can be shown comparing figs. 1 and 2 with fig. 4: 
the EoS parameter presents some oscillations, specially 
for negative redshifts around the phantom barrier, while 
the evolution of the scalar field also oscillates at negative 
redshifts. 




It is well known that f(R) gravity is equivalent to a 
kind of Brans-Dicke theory with a null kinetic term, and 
a scalar potential, (see for example, [7] and references 
therein), 



S = jd 4 x^/^[<j) R-V{<j))+2K 2 Xn] 



(8) 



where the scalar field and the potential are related with a 
particular f(R) action by, 

*=/«(*), V&{R))=f R (R)R-f{R). (9) 

Let's reconstruct the corresponding scalar-tensor theory 
for the particular model considered in the previous sec- 
tion. By the expressions (9), the corresponding relation 
0(7?) for the model (6) is easily obtained, 
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FIGURE 4. Evolution of the scalar field <j>(z), (10). As 
shown, there is a correlation between the evolution of the EoS 
parameter in Figs. 1-2, and the evolution of the scalar field. 

Let us now analyze the phase space {H,<p} for this 
model in vacuum. For simplicity, we assume vacuum 
FLRW equations, so that by combining FLRW equations 
(in the scalar-tensor equivalence), the following equa- 
tions are obtained, 



H = -2H Z + -V' ((b) . 6 = — 
6 yvj 1 v 3H 
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which yield, 



dH 
d(j> 



V't 



2H 2 



(12) 



(13) 



This equation describes the phase space for a particu- 
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FIGURE 5. Phase space H(tj>) for the HS model, where 
V+ in (11) is considered. We have assumed initial conditions 
close to the boundary of the scalar field: (a) 7f(0.95) = H^, 
(b) #(0.95) = 10i?g, (c) H(0.95) = lOQflg. 



lar scalar potential. Fig. 5 the phase space for the model 
(6) with the same free parameters assumed above. Then, 
it shows that regardless of the choice of initial condi- 
tions, the behavior of the Hubble parameter tends to a 
stable point, close to zero, where some tiny fluctuations 
are observed. Other viable models have a similar behav- 
ior, where tiny differences between the models are ob- 
served when the future evolution of the EoS is studied, 
as pointed out in Ref. [5]. 

DISCUSSIONS 

Here, we have analyzed the future evolution of a par- 
ticular viable f(R) gravity model, [3], The cosmologi- 
cal evolution of the effective EoS enters into a phantom 
stage in the present or near future, while the cosmologi- 
cal evolution, similarly to ACDM model at z = 0, is well 
reproduced for small redshifts, where we have imposed 
some particular conditions in order to fit the model with 
the observational constraints. On the other hand, the tran- 
sition to the phantom phase does not imply directly the 
occurrence of a future singularity, since the Hubble pa- 
rameter and its first derivatives remain finite. In addition, 
the scalar-tensor counterpart of the f(R) model shows a 
well behavior of the phase space, at least when (j) < 1, 
which is the case for small redshifts as shown in Fig. 4. 
Nevertheless, this does not prevent that a cosmological 
singularity, induced by the pole of the scalar potential, 
may occur in the past. However, for large redshifts, one 
expects that modifications may be negligible, so that the 
singularity may be avoided. A deeper study on this issue 
should be done in future works. Moreover, the possibility 



of a Little Rip, a phase where the strength of the expan- 
sion would be comparable with some bounded systems 
as the Solar system (see Ref. [9]), seems to be far to occur 
in viable f(R) gravity, as pointed out in Ref. [5], Hence, 
viable modified gravities have still a strong support to be 
considered as a serious candidate to dark energy. 

ACKNOWLEDGMENTS 

I acknowledge support from a postdoctoral contract from 
the University of the Basque Country (UPV/EHU) un- 
der the program "Specialization of research staff", and 
support from the research project FIS201 0-1 5640, and 
also by the Basque Government through the special re- 
search action KATEA and UPV/EHU under program 
UFI 11/55. 



REFERENCES 

1. S. Nojiri and S. D. Odintsov, eConf C0602061, 06 (2006) 
[Int. J. Geom. Meth. Mod. Phys. 4, 115 (2007)]; hep- 
th/0601213; arXiv: 0807.0685; K. Bamba, S. Capozziello, 
S. Nojiri and S. D. Odintsov, arXiv: 1205.3421 [gr-qc]; 

S. Capozziello and M. De Laurentis, Phys. Rept. 509, 167 

(2011) [arXiv: 1108.6266 [gr-qc]]; S. Capozziello and V. 
Faraoni, Beyond Einstein Gravity, Fundamental Theories 
of Physics Vol. 170, Springer Ed., Dordrecht (201 1); A. de 
la Cruz-Dombriz and D. Saez-Gomez, Entropy 14, 1717 

(2012) [arXiv: 1207.2663 [gr-qc]]. 

2. S. Nojiri, S. D. Odintsov and D. Saez-Gomez, Phys. Lett. 
B 681, 74 (2009), arXiv:0908. 1269 [hep-th]. 

3. W. Hu and I. Sawicki, Phys. Rev. D 76 064004 (2007), 
arXiv:0705.1 158[astro-ph]. 

4. S. Nojiri and S.D. Odintsov, Phys. Rev. D 77 026007 
(2008), arXiv:0710.1738[hep-th]; S.A. Appleby and 
R.A. Battye, arXiv:0705.3199[astro-ph]; L. Pogosian 
and A. Silvestri, Phys. Rev. D 77 023503 (2008), 
arXiv:0709.0296[astro-ph]; S. Tsujikawa, Phys. Rev. D 77 
023507 (2008), arXiv:0709.1391[astro-ph]; A. Starobinsky, 
lETPLett. 86 157 (2007), arXiv:0706.2041 [astro-ph]; S. 
Capozziello and S. Tsujikawa, Phys. Rev. D 77 107501 
(2008) arXiv:0712.2268[gr-qc]. 

5. D. Saez-Gomez, arXiv: 1207.5472 [gr-qc]. 

6. E. Komatsu et al. [WMAP Collaboration], Astrophys. J. 
Suppl. 192, 18 (2011) arXiv: 1001.4538 [astro-ph.CO]; 

7. D. Saez-Gomez, Gen. Rel. Grav. 41, 1527 (2009) 
arXiv:0809.1311 [hep-th]; 

8. S. A. Appleby, R. A. Battye and A. A. Starobinsky, ICAP 
1006, 005 (2010) [arXiv:0909.1737 [astro-ph.CO]]. 

9. P. H. Frampton, K. J. Ludwick and R. J. Scherrer, Phys. 
Rev. D 84, 063003 (2011) arXiv: 1106.4996 [astro-ph.CO]; 
S. Nojiri, S. D. Odintsov and D. Saez-Gomez, AIP Conf. 
Proc. 1458, 207 (2011), arXiv: 1108.0767 [hep-th]; 



